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We analyze the entanglement dynamics of a system composed by a pair of neutral two-level atoms 
that are initially entangled, and the electromagnetic field, initially in the vacuum state, within the 
formalism of perturbative quantum field theory up to the second order. We show that entangle- 
ment sudden death and revival can occur while the atoms remain spacelike-separated and therefore 
cannot be related with photon exchange between the atoms. We interpret these phenomena as the 
consequence of a transfer of atom-atom entanglement to atom-field entanglement and viceversa. We 
also consider the different bi-partitions of the system, finding similar relationships between their 
entanglement evolutions. 
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I. INTRODUCTION 

Entanglement between qubits may disappear in a fi- 
nite time when the qubits interact with a reservoir. This 
is commonly known as "entanglement sudden death" 
(ESD) . After its discovery P, 0, y] , the phenomenon has 
attracted great attention (for instance, 3, H, H, H, H, 
[13, El, [13) a.nd has been observed experimentally [l3l |. 

ESD shows up in a variety of systems that can be 
roughly divided in two sets: those in which the qubits in- 
teract individually with different reservoirs and those in 
which they interact with a common environment. In par- 
ticular, in 0, [31 a system of a pair of two- level atoms 
interacting with a common electromagnetic vacuum is 
considered. The dynamics of the system is given in all the 
cases by the Lehmberg-Agarwal master equation [H, [3 
which is derived with the rotating wave approximation 
(RWA) and the Born-Markov approximation. Recently, 
non-Markovian [8] and non-RWA effects have been 
considered in systems of qubits coupled individually to 
different reservoirs. There are good reasons for going be- 
yond the Markovian and RWA scenario in the case of a 
pair of two-level atoms in the electromagnetic vacuum. 
For short enough times non-RWA contributions are rele- 
vant [31 and a proper analysis of causality issues can only 
be performed if they are taken into account [3, [13, IMf ■ 
Besides, as we shall show in this paper the death of the 
entanglement between the atoms is related with the birth 
of entanglement between the atoms and the field, and 
therefore the field is actually a non- Markovian reservoir. 
This was also the case in [3, [3, [HI with different reser- 
voirs. 
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turbative quantum electrodynamics (QED) to the sys- 
tem of a pair of neutral two-level atoms interacting lo- 
cally with the electromagnetic field, and for initially sep- 
arable states analyzed the generation of entanglement. 
This is a non-Markovian, non-RWA approach. The use 
of the Lehmberg-Agarwal master equation can be seen as 
a coarse-grained in time approximation to the perturba- 
tive treatment ^24] . The first goal of this paper is to apply 
also the QED formalism to analyze the ESD in these sys- 
tems for initially entangled atomic states^omparing the 
results with the previously obtained i\ with master 
equations. We will focus mainly on the range r/{ct) « 1, 
r being the interatomic distance and t the interaction 
time, in order to investigate the role of locality. We will 
also consider for the first time in these systems the rest of 
pairwise concurrences, namely the entanglement of each 
atom with the field, and multipartite entanglement, fol- 
lowing the spirit of [3, 3, [3, [3| • While the mentioned 
papers deal with a four qubit model, our model here con- 
sists in two qubits (the atoms) and a qutrit (the electro- 
magnetic field, which may have 0, 1 or 2 photons). We 
shall show that the phenomenon of revival of entangle- 
ment after the ESD [7,] can occur for r > ct, and therefore 
is not related with photon exchange as is usually believed. 
We will see that atom-atom disentaglement is connected 
with the growth of atom-field entanglement and vicev- 
ersa. Similar relationship will be obtained among the 
"atom-(atom-l-field)" and "field- (atom-f atom)" entangle- 
ments. 

The reminder of the paper is organized as follows. In 
section II we will describe the Hamiltonian and the time 
evolution from the initial state of the system. In section 
III we will obtain the reduced state of the atoms and 
analyze the behavior of its entanglement. In section IV 
the same will be performed with the reduced state of each 
atom and the field, comparing the entanglement cycle 
with the one obtained in the previous section. Tripartite 
entanglement will be considered in section V in terms of 
the entanglement of all the different bi-partitions of the 
system, and we conclude in section VI with a summary 
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of our results. 



II. HAMILTONIAN AND STATE EVOLUTION 

To address the atom- field interactions, we assume that 
the relevant wavelengths and the interatomic separation 
are much larger than the atomic dimensions. The dipole 
approximation, appropriate to these conditions, permits 
the splitting of the system Hamiltonian into two parts 
H = Ho + Hi that are separately gauge invariant. The 
first part is the Hamiltonian in the absence of interac- 
tions other than the potentials that keep A and B sta- 
ble, Hq = Ha + Hb + -ff field- '^^^ second contains all 
the interaction of the atoms with the field, which in the 
dipole approximation we will use is given by: 

Hj = -- y d„(x„,i)D(x„,i), (1) 

where D is the electric displacement field, and d„ = 
J2i e / d^Xi {E\{xi — x„) I G ) is the electric dipole mo- 
ment of atom n, that we will take of equal magnitude 
for both atoms (d dA = de) as required by angular 
momentum conservation in the photon exchange. | E ) 
and I G ) will denote the excited and ground states of the 
atoms, respectively. 

In what follows we choose a system given initially by 
an atomic entangled state, with the field in the vacuum 
state I 0) : 

\i,)o ={a\EE) + (3\GG))-\0). (2) 

The system then evolves under the effect of the interac- 
tion during a lapse of time t into a state: 

|^)t=r(e-^/o^*'^^(*')/^)|V)o, (3) 

T being the time ordering operator. Up to second order 
in perturbation theory, ([3]) can be given in the interaction 
picture as 

|atoml,atom2, field)* = a\EEO)t + /3|GG0)t (4) 
where 

\EEO)t= {{l + a)\EE) +b\GG))\0) 
+ {uA\GE)+UB\EG))\l) + {f\EE)+g\GG))\2){5) 

and 



where 

b = (0|T(55 5;,)|0),6' = (0|T(5+5+)|0), 

UA = (l|5^|0),i;^ - (1|5+|0) (7) 

UB - {l\Ss\0),VB = (1|5+|0) 

f = l9{h-t2){2\S^{h)S+{t2)+Ss{h)S+it2)m, 
9 = {2\T{S-sS-^m, g' = {2\T{S+S\m 

being S = -fj^ dtHi{t') = S+ + S~ , T the time or- 
dering operator and \n), n = 0, 1, 2 is a shorthand for 
the state of n photons with definite momenta and polar- 
izations, i.e. 1 1) = |k, e), etc. Here, a and a' describe 
intra-atomic radiative corrections, ua{ub) and va {vb) 
single photon emission by atom A {B), and g and g' by 
both atoms, while / and /' correspond to two photon 
emission by a single atom. Only b and b' correspond to 
interaction between both atoms. The sign of the super- 
scripts is associated to the energy difference between the 
initial and final atomic states of each emission or absorp- 
tion. In Quantum Optics, virtual terms like a', &, b' , va, 
''^B, fi f and g' , which do not conserve energy and ap- 
pear only at very short times, are usually neglected by 
the introduction of a RWA. In the dipole approximation 
the actions hS^ in ([5]) reduce to 

^ f dt'e±'"*'dE(x,t') (8) 
ft Jo 

where £7 = uje — is the transition frequency, and we 
are neglecting atomic recoil. This depends on the atomic 
properties 17 and d, and on the interaction time t. In 
our calculations we will take (17|d|/ec) = 5 • 10^"^, which 
is of the same order as the Is — > 2p transition in the 
hydrogen atom, consider S7i > 1, and focus mainly on 
the cases (r/ct) ~ 1. Therefore, | i?) is actually a triply 
degenerate state \E,m) with m = 0,±1 and we will 
average over two different independent possibilities for 
dipole orientations: d^ = ds = d = du^ for transitions 
with Am = [25] and d = d{ux ±iUy)/V2 [25,] for 
transitions with Am = ±1. 



III. SUDDEN DEATH AND REVIVAL OF 
ATOM-ATOM ENTANGLEMENT 



1 GGO)f - ((1 + a') 1 GG) -I- b' \EE))\ 0) After tracing over all the states of the field, the density 

+ {va\EG) + vb\G E))\l) + {f \GG) + g' \E E))\2){6) matrix of the atomic state pab takes the form (in the 
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basis {\EE), \EG), \GE), \GG)}): 



and 



/ pii pii 

P22 P23 

P*23 P33 

V P44 



(9) 



where 



pii = |a(l + a) + + \af + ^ g' \\ 

P22 = P33 - + |/3|2|«|2+2i?e(a/3*r) 

P44 = |a6 + + a')l' + l«5 + (10) 

P14 = |a|2((l + a)&* +/.9*) + |/3|2((l + a')*fc' 

+ g' n + a(3*{(l+a){l + a') + / /'*) 

P23 = \a\'^UBU*A + \(3\'^VAV*B + 2Re{aP*uv*) 

N = Pll + P22 + P33 + PAi 

2 _ I 12 I 12 _ U, , 12 _ U,„ 12 



where | w p 
I 



= ~ ubv\ and w v* = u\ v\ = ub Vg. 

The computation of a, &, etc. can be performed fol- 
lowing the lines given in the Appendix A of [2^ . where 
they were computed for the initial state \E G) and only 
for Am — 0. In terms of z = ^r/c and x = r/cA, being 
r the interatomic distance, we find: 



(ln(l - 
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with K = a\d |V(e2 r^) and / = /+ 



in), 

(11) 

/_ , where: 



I± 



2z 



±2co8(-) e^'' Ei(Tiz)+ e"''^^^i'^ 

X 

■■(^^i^ Ei{-iz{l±-))]{12) 



Ei(iz(l± -)) - e' 

X 



X 

iz (l-l/x) 



for a; > 1, having the additional term —2 Trie 
otherwise. 

jwp, Ijiif, I, ubu\, vaVb and uv* have been com- 
puted in [22] ■ Besides: 

g = ubu'^ + uau'b , g' = vav'b + vb v'a 

f = e{ti-t2){vA{tl)u'A{t2) + UA{tl)v'A{t2) 

+ VB{h)u'B{t2) + UB{h)v'B{t2)) (13) 

/' = e{tl-t2){uA{tl)v'j^{t2) + u'j,{ti)vA{t2) 

+ UB{ti)v'B{t2) + u'B{ti)vB{t2)) 

where the primes are introduced to discriminate between 
the two single photons. 

We will use the concurrence C(/9) |26j| to compute the 
entanglement, which for a state like © is given by, if 

|P14 I 



VP22 P33 + 
'C{pAB 



P23 I > VPiI P44 

■2(|p23| - VPll P44) 







(14) 



C(/9ab) = max 



2(|pi4| - VP22P33) 



(15) 



otherwise. 



If we take a = ^Jp and [3 = yT 
appears at a range of values of 
creasing r, in agreement with 



~p, we find that ESD 
that decreases with in- 
Although this would 
suggest that ESD disappear for r large enough, we find 
that there are high values of p for which ESD exists for 
arbitrary large r. In Fig. 1, we represent C{pab) in 
front of X for different values of z and p = 0.98. ESD 
occurs at z/x = fit of the order of 10''. Thus, as z (that 
is r) grows, ESD is shifted to higher values of x. It is 
also interesting to analyze the phenomenon of entangle- 
ment revival, discovered in these systems in [7]. We find 
that the dark periods between death and revival has 
larger time durations for increasing z. Besides, although 
in [TI the revival is described as a consequence of the pho- 
ton exchange, for r sufficiently large both the ESD and 
the revival can occur for a; > 1, where photon exchange 
is not allowed. We think that the explanation for en- 
tanglement revival is closer to the spirit of 0, [ll| where 
entanglement revival between noninteracting atoms is in- 
terpreted as coming from entanglement transfer between 
different parts of the system. We shall discuss this point 
in the following sections. 
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FIG. 1: Concurrence C(pab) in front of a; = r/ct for p = 0.98 
and z = Slr/c = 2 ■ lO'' (solid line), 5 • 10*^ (dashed line) and 
2 • 10^ (dotted line). In the latter case sudden death and 
revival of entanglement occur for a; > 1. 



In Fig. 2 we sketch the dependence with p. Although 
sudden death and revivals appear in a very restricted 
range of the parameter, they are only a particular case of 
the generic behavior of entanglement observed in a wider 
range, which can be described as disentanglement up to 
a minimum value and growth of quantum correlations 
since then. 
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FIG. 2: Concurrence C(pAs) in front of a; = rjct for z = 
VLrjc = 2-10^ and p = 0.97 (solid line), p = 0.98 (dashed line) 
and p = 0.99 (dotted line). In the first case, entanglement 
decreases as t grows up to a minimum value and begin to 
grow since then. This behavior becomes entanglement sudden 
death and revival when the minimum value is for higher 
values of p. C{pab) tends to as s — + cx) and p — + 1. 



IV. ATOM-FIELD ENTANGLEMENT 

Tracing Q over states of atom A (B) the reduced 
atom- field density matrix pbf (paf) is obtained. Tak- 
ing the basis {| £; ), | 1 ), | £; 2 ), | GO ), | G 1 ), | G 2 )}, 
we have: 



Pbf = pAF 
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(16) 



with 

p'll = 

P33 = 

P'm = 

Pis = 

P24 = 

P'26 = 

P35 = 

P46 = 

iV' = 



partial transposes of both pbf and p^j? are 

p'll + P55 ± V(p'll-P55)'+4|pi4P 



|a(l + a) + /3 5'|2, p'^^ = p'^^ = 
l«/ + /35'l',P;4 = l"6 + /3(l + a')l' 
I « .9 + /3 /' P'i3 = (« (1 + «) + /3 &') (a / + /35') 
(a(l + a) -H/3fe') (qub + I3va)* 
{auA + I3vb) iP (1 + a') + ab)* 
{auA + f3vB){ag + (3 f')* 
{af + fig') (auB + I3va)* 
[ab + f3{l +a'))iag + /?/')* 



(17) 
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P44 



P'55 



P66 



There are no operational generalizations of concurrence 
for mixed states in 2 x 3 dimensions like the ones in Eq. 
(fTH]) . We will use the negativity [27i| N(p), which is the 
absolute value of the sum of the negative eigenvalues of 
the partial transposes of a state p. For the 2x2 and 2x3 
cases N(p) > is a necessary and sufficient condition for 
p to be entangled. 

Up to second order in perturbation theory, we have 
that N' = N and that the nonzero eigenvalues of the 



2N' 



and 



P44 + P55 ± ViP'44 - P22)' + 4|pi5p 
2N' 



(18) 



(19) 



being zero the other two. In Eqs. and (fTO|) only 

the terms up to second order are retained. Therefore, if 
IP24P > p'll P55 tJ^en A_ < and if \p[^\^ > p'22 P44 then 
A'_ < 0. 

In Fig. 3 we represent N{pbf) = ^(paf) in front of x 
for same values of p and z of Fig. 1. We see that the neg- 
ativity grows from at a: — > c» (t = 0) to its maximum 
value and then starts to decrease and eventually vanishes, 
following the opposite cycle to the entanglement of pab ■ 
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FIG. 3: Negativity N{pbf 

p = 0.98 and z ^ nr/c ^ 2 ■ 10** (solid line), 5 ■ 10'' (dashed 
line) and 2 • 10^ (dotted line). Entanglement increases from 
at a: ^ 00 up to a maximum value and then decreases and 
vanishes eventually. 

Although it would be interesting to look for conser- 
vation rules of entanglement like the ones in |Tl|, [l^] , 
this search is beyond the focus of this paper since in our 
study we are using different entanglement measures in 
Hilbert spaces of different dimensions. Besides, except 
for the concurrence between atoms A and B, the rest of 
the concurrences in the mentioned papers have not ob- 
vious counterparts in our case. But it is clear that in 
general the entanglement cycle between atoms is corre- 
lated with the entanglement cycle between each atom and 
the field, as can be seen in Fig. 4 in a particular case. 
Although atom-field entanglement may change while the 
other remains zero, both entanglements cannot increase 
or decrease at the same time. 



V. TRIPARTITE ENTANGLEMENT 

Tripartite entanglement has been widely studied in 
terms of the entanglement of the different bipartitions 
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FIG. 4: Negativity N{pbf) ~ N(paf) (solid line) and concur- 
rence C{pab) (dashed line) in front of a; = r/ct for p — 0.98 
and z = Qr/c = 5 ■ 10®. Entanglement atom-atom cycle 
is clearly correlated with the atom-field cycle, although the 
sum is not a conserved quantity. Although atom-field en- 
tanglement may change while the other remains zero, both 
entanglements cannot increase or decrease at the same time. 



A - BC, B - AC, C ~ AB in the system |28|, \M l30|, 
where A, B and C stand for the three parties. Here, 
we will compute the I concurrences [3l| 'Ca-bf, Cb_af, 
Cf-ab, where Cj^kl = ■\/2(l — Tr p'j), where J runs 
form ^ to F and KL from BF to AB respectively, being 
Pj the reduced density matrix of J. A and B stand for 
the atoms, and F for the field. 
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FIG. 5: I concurrence Ca-bf ~ Cb-af in front of a; = r/ct 
for p = 0.98 and z = nr/c = 2 ■ 10^ (solid line), 5 • 10^ (dashed 
line) and 1 ■ 10® (dotted line) . Entanglement disappears faster 
than the entanglement between the atoms (Fig.l) and remains 
since then. 

Tracing ^ over BF (AF), we find the following den- 
sity matrices pA (pb)- 



PA = Pb 



1 



PAll 
PA22 



(20) 



where pAii = p'n + P33 + P22 and pA22 = p'm + Pee + 
P22 and Na = Pah + Pa22- In Fig- 5 we sketch the 
behavior of Ca-bf and C-b-af m front of x for different 



I concurrence Ca-bf ~ 

l6 



-AF in front of a:: = 
r/ct for z — Q,r/c = 2 ■ 10® and p = 0.50 (solid line), p — 
0.75 (dashed line) and p = 0.98 (dotted line). Entanglement 
sudden death occurs for a wider range than the entanglement 
between the atoms (Fig. 2). 



values of z. Entanglement vanishes before the death of 
the entanglement between A and B, and does not have 
a revival. Besides, ESD appears in a wider range of p, as 
can be seen in Fig. 6. 
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FIG. 7: I concurrence Cf-ab in front oi x — r/ct for p = 0.98 
and z = fir /c = 2 ■ 10^ (solid line), 5 ■ 10^ (dashed line) 
and 1 • 10® (dotted line). Entanglement grows from to its 
maximum value at a; « 0.1 and then decreases. 



Now, tracing ^ over AB we obtain the reduced den- 
sity matrix of the field pp- 



PF 




Pfh Pfi3 

PF22 

PF33 



(21) 



where pFii = p'n + p'44, PF22 = 2p22, PF33 = p'33 + Pea, 
PF02 = p'i3 + Pie and Np = PFii + PF22 + P_F33- In 
Fig. 7 we represent Cf-ab m front of x for the same 
values of z and p as in Fig. 5. Entanglement grows from 
to a maximum value at a; ~ 0.1 and then decreases. 
The growth of Cf~ab is correlated with the decrease of 
Ca-bf and Cb-af in the same way as the magnitudes 
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FIG. 8: I concurrence Cf-ab (solid line) and Ca-bf = 
Cb-af (dashed line) in front of a; = r/ct for z = Qr/c = 
2 • 10^ and p = 0.98. Both magnitudes cannot increase or 
decrease at the same time. 



analyzed in the previous section, as can be seen in Fig. 
8 for a particular case. 



sisting in a pair of neutral two-level atoms A and B inter- 
acting with a common electromagnetic field F. At t = 
atoms are in the Bell state \E E) + ^1 — p \GG) 
and the field in the vacuum state. The evolution of 
this state has been considered within the non-Mar kovian, 
non-RWA approach of quantum electrodynamics up to 
second order in perturbation theory. We find ESD and 
revival of entanglement in the reduced state of the atoms, 
in a range of p that decreases with the interatomic dis- 
tance r, in agreement with the results obtained with mas- 
ter equations 0. For r large enough, we find that the 
revival of entanglement can occur with r > ct and there- 
fore is not a consequence of photon exchange between the 
atoms. We find that this phenomenon is strongly related 
to the transfer of entanglement between the different sub- 
systems of two parties that coexist in the entire system: 
we obtain sort of entanglement cycle for the atom-field 
reduced states opposite to the atom-atom one. We have 
considered also the different bi-partitions of the system, 
namely A ~ BF, B — AF and F — AB, finding similar 
relationships between their entanglement cycles. 



VI. CONCLUSIONS 



Acknowledgments 



We have analyzed in a previously unexplored space- 
time region the entanglement dynamics of a system con- 



This work was supported by Spanish MFC FIS2005- 
05304 and CSIC 2004 5 OE 271 projects. 



[1] K. Zyczkowski, P. Horodecki, M. Horodecki, and 

R. Horodecki, Phys. Rev. A 65, 012101 (2001). [18 
[2] L. Diosi, Lect. Notes Phys. 622, 15 (2003). [19 
[3] T. Yu and J. H. Eberly, Phys. Rev. Lett. 93, 140404 

(2004). [20 
[4] T. Yu and J. H. Eberly, Phys. Rev. Lett. 97, 140403 

(2006). [21 
[5] M. Yonag, T. Yu, and J. H. Eberly, J. Phys. B: At. Mol. 

Opt. Phys. 40, S45 (2007). [22 
[6] A. Jamroz, J. Phys. A: Math. Gen. 39, 7727 (2006). [23 
[7] Z. Ficek and R. Tanas, Phys. Rev. A 74, 024304 (2006). [24 
[8] B. Bellomo, R. L. Franco, and G. Compagno, Phys. Rev. 

A 77, 032342 (2008). 
[9] J. P. Paz and A. J. Roncagha, Phys. Rev. Lett. 100, [25 

220401 (2008). 

[10] C. E. Lopez, G. Romero, F. Lastra, E. Solano, and J. C. [26 

Retamal, Phys. Rev. Lett. 101, 080503 (2008). 
[11] S. Chan, M. D. Reid, and Z. Ficek, arXiv[quant- [27 

ph]:0810.3050. 

[12] J. H. Cole, arXiv[quant-ph]:0809.1764. [28 
[13] M. P. Almeida, F. de Melo, M. Hor-Meyll, A. Salles, S. P. [29 
Walborn, P. H. S. Ribeiro, and L. Davidovich, Science 
316, 579 (2007). [30 
[14] Z. Ficek and R. Tanas, Phys. Rev. A 77, 054301 (2008). 
[15] R. H. Lehmberg, Phys. Rev. A 2, 883 (1970). [31 
[16] G. S. Agarwal, Springer Tracts Mod. Phys 70, 1 (1974). 
[17] F.-Q. Wang, Z.-M. Zhang, and R.-S. Liang, arXiv[quant- 



ph]:0806.3306. 

J. Leon and C. Sabm, Phys. Rev. A 78, 052314 (2008). 
E. A. Power and T. Thirunamachandran, Phys. Rev. A 
56, 3395 (1997). 

P. W. Milonni, D. F. V. James, and H. Fearn, Phys. Rev. 
A 52, 1525 (1995). 

A. K. Biswas, G. Compagno, G. M. Palma, R. Passante, 

and F. Persico, Phys. Rev. A 42, 4291 (1990). 

J. Leon and C. Sabm, arXiv[quant-ph]:0804.4641. 

J. Leon and C. Sabm, arXiv[quant-ph]:0805.2110. 

C. Cohen-Tannoudji, J. Dupont-Roc, and G. Gryn- 

berg. Atom-photon interactions (Wiley Interscience, New 

York, 1998). 

P. W. Milonni and P. L. Knight, Phys. Rev. A 10, 1096 
(1974). 

S. HiU and W. K. Wootters, Phys. Rev. Lett. 78, 5022 
(1997). 

G. Vidal and R. F. Werner, Phys. Rev. A 65, 032314 
(2002). 

H. -S. Song and C.-s. Yu, Phys. Lett. A 330, 377 (2004). 
C. Sabm and G. Garci'a-Alcaine , Eur. Phys. J. D 48, 
435 (2008). 

P. Facchi, G. Florio, and S. Pascazio, Phys. Rev. A 74, 
042331 (2006). 

P. Rungta, V. Buzek, C. M. Caves, M. Hillery, and G. J. 
Milburn, Phys. Rev. A 64, 042315 (2001). 



